Magnetic-Field-Induced 4/-Octupole in CeBe Probed by Resonant X-ray Diffraction 
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CeBe, a typical Fg-quartet system, exhibits a mysterious antiferroquadrupolar ordered phase in 
magnetic fields, which is considered as originating from the T^^y^-type magnetic octupole moment 
induced by the field. By resonant x-ray diffraction in magnetic fields, we have verified that the 
Txyz-type octupole is indeed induced in the 4/-orbital of Ce with a propagation vector (|, |, 
thereby supporting the theory. We observed an asymmetric field dependence of the intensity for an 
electric quadrupole (E2) resonance when the field was reversed, and extracted a field dependence 
of the octupole by utilizing the interference with an electric dipole (El) resonance. The result is in 
good agreement with that of the NMR-line splitting, which reflects the transferred hyperfine field at 
the Boron nucleus from the anisotropic spin distribution of Ce with an O^^y-type quadrupole. The 
field-reversal method used in the present study opens up the possibility of being widely applied to 
other multipole ordering systems such as Np02, Ce^^Lai-j^Be, SmRu4Pi2, and so on. 

PACS numbers: 61.05.cp , 71.27.-|-a , 75.25.+Z , 75.10.Dg 



A rich variety of electronic phases arising from mul- 
tiple degrees of freedom of / electrons has attracted 
great interest in recent years. In addition to mag- 
netic dipole moment, electric quadrupole and magnetic 
octupole, etc., behave as independent degrees of free- 
dom in a crystal-field eigenstate with orbital degener- 
acy. Quadrupole orders are frequently realized in local- 
ized /-electron systems and, more exotically, octupole 
orders can also takeplace as in Np02 and Cex^s.i-xQ& 
[x < 0.8) [l|,Sii,i. Furthermore, it has recently been 
recognized that these multipoles sometimes play fun- 
dametal roles when / electrons are hybridized with con- 
duction electrons. In Pr-based filled skutterudites such 
as PrRu4Pi2, a 4/-hexadecapole order is combined with 
a Fermi-surface nesting, causing a metal-insulator transi- 
tion @. In PrOs4Sbi2, it is suggested that a quadupolar 
excitation is associated with the heavy-fermion supercon- 
ductivity Thus, understanding the physics of multi- 
pole moments is of fundamental importance. 

One of the difficulties of multipole research is that they 
are hard to identify as is often expressed as hidden or- 
der paremeter. In most cases, a primary order parameter 
(OP) is initially inferred indirectly by combining various 
pieces of information from macroscopic and microscopic 
methods. Then, detailed investigation of a secondary OP 
by neutron and x-ray diffraction, and of a hyperfine field 
by NMR, using a single crystal, may provide evidence 
for the multipole OP [l|, Among these microscopic 
probes, resonant x-ray diffraction (RXD) has a distinc- 
tive ability to directly probe ordered structures of mul- 
tipole tensors up to rank 4, using an electric quadrupole 
(E2) resonant process Q. With respect to the observa- 
tion of an antiferroquadrupole (AFQ) order, there have 
already been several examples of successful applications 
of RXD mainly using^an electric dipole (El) resonance, 
typically for DyB2C2 [9|. On the other hand, with respect 
to an antiferrooctupole (AFO) order, there has been only 



one report on Ceo.7Lao.3B6 by Mannix et al. [4]. They 
successfully detected an E2 signal at zero field, mea- 
sured the azimuthal-angle dependence, and concluded an 
AFO order, which was also confirmed by neutron diffrac- 
tion 0. However, it was pointed out that the contribu- 
tion from the 4/-quadrupole to the E2 signal cannot be 
ruled out [l3|. Azimuthal dependence only is not suffi- 
cient to separate contributions from different rank tensors 
to an E2 signal. 

In this Letter, we report an effective method that can 
distinguish between even and odd rank tensors, which 
will be quite useful in studying octupole orders, especially 
those induced in magnetic fields. Since various kinds of 
multipoles are induced in magnetic fields and affects the 
macroscopic properties, it is of fundamental importance 
to trace what kind of multipole is induced in magnetic 
fields. 

A compound we study is CeBg, a typical Fg-quartet 
system with a simple cubic structure. The Fg has 15 
degrees of freedom in total, 3 dipoles, 5 quadrupoles, 
and 7 octupoles ll|. At zero field, CeBg exhibits an 
Oa;y-type AFQ order at rQ=3.3 K followed by an an- 
tiferromagnetic (AFM) order at Tn=2.3 K [H, [H, [H | . 
In magnetic fields, Tq exhibits an anomalous increase 
up to 8.3 K at 15 T [ij], whose most important mech- 
anism has been ascribed to an antiferro-type interaction 
between field-induced octupoles of T^yz-^ypQ [HI- Split- 
ting of the Boron-NMR line in the AFQ phase can be 
a strong evidence for this interpretation [3, It is 
explained by a phenomenological analysis of the hyper- 
fine field at the Boron nucleus in terms of the multi- 
pole moments of Ce based on symmetry arguments. To 
be exact, however, we have to mention that direct ev- 
idence for the existence of octupole is still lacking. As 
pointed out by Hanzawa, the microscopic mechanism of 
the NMR splitting is due to the transferred hyperfine field 
(THF) via the 2p and 2s conduction electrons, reflect- 
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FIG. 1: top: Fluorescence spectrum of CeBe. bottom: X- 
ray energy spectra of the (| | ^) superlattice reflection in 
magnetic fields with reversed directions. The triangle shows 
the background due to the fluorescence. 

ing the anisotropic spin distribution of Ce with an O^y- 
type quadrupole [17] . The splitting can also be explained 
phenomenologically if one considers that the THF, e.g., 
for Borons along the z-axis, is proportional to mF(Oxy) 
where mp is a field-induced uniform magnetization [l8|. 
However, mp{Oxy) is not identical to Txyz in the sense 

that Txyz = ^{JxOyz+JyOzx+JzOxy) rcprescnts a com- 
plex magnetization distribution where all the three terms 
are equally induced even for the field along the z-axis. 
The field-induced octupole is essentially a quantum me- 
chanical phenomenon and the increase of Tq in magnetic 
fields requires consideration of this real Txyz 4/-octupole. 
Therefore, it is worth verifying whether the Txyz-octwpo\e 
is induced in the Ce 4/-orbital itself. 

RXD experiment has been performed at Beamline 3A 
of the Photon Factory in KEK, using a vertical field su- 
perconducting magnet on a two-axis diffractometer. A 
sample with a mirror-polished (331) surface was mounted 
in the cryostat, so that the [001] and [110] axes were in 
the horizontal scattering plane and the field was along 
the [110] axis. The incident photon was 7r-polarized and 
the energy was tuned to the Ce Lni absorption edge. 

Figure [1] shows the energy spectra of the (| | i) su- 
perlattice reflection at 2.5 K in the AFQ phase for several 
magnetic fields with reversed directions. We immediately 
notice that the peaks at 5.724 keV (El) and at 5.718 keV 
(E2) become stronger and well resolved for fields in the 
plus direction, whereas for fields in the minus direction 
the E2 peak becomes obscure. From this result, we can 



extract the field dependence of the quadrupole and oc- 
tupole moments as explained next. 

The energy- and field-dependent structure factor for 
resonant diffraction is generally expressed as 

F,.,so{E,H) =ZEiiH){fUE)+ifi^AE)} 

+Zp2{H){f{,^{E)+if^^{E)}, (1) 

where Zei and Ze2 are unit-cell structure factors for 
El and E2 processes, which are directly coupled with 
the atomic tensors (Tg^^) of hd and 4/ orbitals, respec- 
tively [8|]. They are written as 

n K=Q q 

Ze2 =5:e-«" ^ ^i?^(rf ))W)<)(-1)«. (3) 

n K=0 q 

Here, Ak and Bk are constant factors for the r&nk-K 
terms, AT^-^) and i?^^) are spherical tensors of the x-ray 
beam determined by the diffraction geometry, i?„ is a 
position vector of the nth Ce ion in a unit cell, and n is 
a scattering vector. {Tq^'^) varies with the applied field. 
The energy dependent term in Eq. ([T]) can be written 
as f{E) = 1/(S - A + ir/2) (A = Aei or Ae2) when a 
resonance can be modeled by a single oscillator. However, 
we leave it here as f '{E)+if"(E) because the actual form 
is not such simple [10| . 

The asymmetry with respect to the field reversal can 
be understood by considering the following two effects. 
The first is that the El and E2 terms interfere, i.e., the 
intensity is proportional to [Zei/ei + ^E2/e21^ and not 
to [Zei/ei]^ + 1^E2/e21^- The second is that the odd 
rank tensor (magnetic dipole and octupole) reverses its 
sign with the field reversal, whereas the even rank ten- 
sor (electric quadrupole and hexadecapole) do not change 
sign. That is, the even(odd) rank terms in Z are sym- 
metric(asymmetric) with respect to the field reversal. To 
analyze the symmetry and asymmetry of the intensity, 
we write the Z factor in Eq. ([T]) as Z^ + iZ'^, where 
Z^{Z^) represents the symmetric(asymmetric) part cor- 
responding to the even(odd) rank term. It is noted 
that the odd rank term is imaginary. The energy and 
field dependent intensity I{E, H) can be calculated by 
\F-[cso{E , H)\'^ , and the symmetric and asymmetric part 
of the intensity, P{E, H) and I'^{E, H), are obtained by 
{I[E, H) + I{E, -H)}/2 and {I{E, H) - I{E, -H)}/2, 
respectively. They are expressed as 

PiE, H) ^{{Zl,,f + {Z^,fm^,f + {f^,f} 

+ {(Z|2)' + (^i2)n{(/E2)' + (/E2)n 

+ 2(Z|iZ|2 + ^ii^g2)Re{/^i/E2} (4) 
nE, H) ^2{Zl^Zl^ - Zl^Zl^l^aiJ^^}. (5) 
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CeBs, (3/2, 3/2, 1/2), T=2.5 K, H//[-1 1 0] 
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FIG. 3: X-ray energy dependence of the symmetric and asym- 
metric components deduced from the data for ±2 T. The solid 
lines are the fits to the data with Lorentzian components. The 
single dotted, double dotted, and dotted lines represent 1st, 
2nd, and 3rd term in Eq. respectively. 



FIG. 2: top: Magnetic-field dependence of the integrated in- 
tensity of the (| I i) resonant Bragg diffraction measured 
at E = 5.724 keV (left) and 5.718 keV (right), corresponding 
to the El and E2 processes, respectively, bottom: Symmetric 
(circles) and asymmetric (squares) components of the inten- 
sity for the respective processes. 



In Fig. [21 we show the field dependence of the inte- 
grated intensity for a rocking scan at each resonance en- 
ergy. The symmetric and asymmetric components de- 
duced from the raw data are shown in the bottom fig- 
ures. At E = Aei the asymmetric intensity is negligibly 
small, whereas it clearly exists at i? = Ae2. From these 
data and Eqs. ^ and ([5]), the field dependence of the 
multipole tensors can be extracted. Of course, to deter- 
mine all the parameters in general, we need information 
from azimuthal-angle dependence, polarization analysis, 
model calculation, and also from other experimental re- 
sults. In the present case of CeBg for || [1 1 0], however, 
some factors can be neglected and the situation become 
quite simple and suited for a demonstration. 

The sharp anomaly in intensity around 0.1 T corre- 
sponds to the one reported in [l9| . This is a phase tran- 
sition from the {Oxy)-OP at zero field, with (Oyz) and 
{Ozx) domains equally populated, to the {aOyz+ f30zx + 
jOxy)-OP, where (a, (3, 7) is the unit vector of the field 
direction. This has also been observed by non-resonant 
x-ray diffraction (20| . Although this is also an important 
nature of the AFQ phase of CeBg, we do not dealt with 
it because it is outside the subject of this Letter. 

Figure [3] shows the energy spectra of P{E) and I^{E) 
deduced from the data for ±2 T. We observe that P{E = 
Aei) is dominated by the first term in Eq. In fitting 
P{E) and P{E), every term in Eqs. ^ and ^ was 



assumed as a Lorentizan, where the real and imaginary 
parts of /ei/e2 connected by the Kramers-Kronig 
relation. Absorption effect and a Gaussian resolution of 
2 eV were also taken into account in the fit. Although 
there are two contributions from Z^^^ and Z^^ to P{E — 
Aei), ^iii' reflecting the fleld- induced AFM dipole, can 
be neglected here. This is justified by the variation of 
I{E ~ Aei) as measured by rotating the crystal around 
the [331] axis. The result can be perfectly explained by 
considering only the AFQ-OP of {aOyz + (30zx + jOxy), 
indicating negligible contribution from the induced AFM. 
In addition, below Tn, we could not detect any signal at 
superlattice spots of the AFM order such as (| | i), 
probably because it was too small. The dipole moment 
in the AFM phase estimated by neutron diffraction is 
^0.28 whereas that for the induced AFM in the 

AFQ phase is -0.05 /ib at H=2 T 18]. That is, the 
dipole is not the main polarization of the 4/ shell giving 
rise to the resonant signal. 

By taking the square root of P{E — Aei), -^ei' 
fleeting the AFQ moment of the 5d orbital, is obtained. 
This is proportional to that of the 4/ orbital, which is 
{Oyz — Ozx) for -ff ]] [1 1 0] from the structure-factor 
calculation. Next, since Z^^ in Eq. ^ can be neglected, 
we can deduce Z^2 by dividing I'^{E = Ae2) by 
As described above, the dipole contribution to Z^^ ™ay 
also be neglected, the obtained result is considered as re- 
flecting only the 4/-octupole. From the structure-factor 
calculation, Z^2 is proportional to {Txyz + 0.02Tf„) for 
H \\ [1 10], where (Txyz) is dominant. These results 
are plotted in Fig. IH ZI2, reflecting the 4/-quadrupole 
and hexadecapole with the same symmetry of Tj'^ — r|g , 
can also be deduced after some data treatments, but this 
results in the same field dependence as that of Z^-^ as 
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the field-reversal method used in the present study can 
be widely applied to other multipole ordering systems 
such as Np02, Cex^a^i-x^e, SmRu4Pi2 and so on. 
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FIG. 4: Magnetic-field dependences of the AFQ and AFO 
moments deduced from the symmetric El and asymmetric 
E2 intensities in Fig. [5] The solid line is a guide for the eye. 
The crosses represent the THF as deduced from NMR. 



expected. 

In Fig. 31 the field dependence of the 4/-octupole shows 
a good agreement with that of THF at the Boron site 
as deduced from NMR 3, 21 1. In addition, it exhibits 
a convex dependence like a Brillouin function. This is 
quite a contrast to the concave field dependence of the 
induced AFM as measured by neutron diffraction [23 |. 
This fact also supports that is dominated by the 
octupole contribution. 

One of the reasons we could obtain this elegant result 
is that the scattering geometry for H \\ [110] provides an 
ideal situation. Using the wave-functions obtained from 
a realistic mean-field model [23| , we can calculate the Z- 
factors for each rank and polarization channel. Firstly, 
all the Z-factors for the tttt' channel vanishes except for 



^E2 tttt' ' ina-king the measurement and analysis straight- 
forward. Vanishing of the signal for the tttt' channel was 
checked by a polarization analysis using a Mo-(200) crys- 
tal analyzer at H=+2 T, though the data in this paper 
were taken without analyzing the polarization. Secondly, 



Z^^^^, and Z^2Tra' ^^i^*^ their maximum at || [1 1 0], 
giving rise to the strongest interference. 

To summarize, we have demonstrated that even and 
odd rank multipoles can be extracted effectively by mea- 
suring the asymmetrical intensity of RXD with respect 
to the field reversal, originating from the interference be- 
tween the El and E2 resonances. In the present case for 
CeBg in H II [1 1 0], this method was quite effective to 
extract the field dependeces of AFQ and AFO moments. 
The result for octupole showed a good agreement with 
that of THF at the Boron site deduced from NMR. Our 
observation directly shows that the octupole moment is, 
indeed, induced in the 4/ orbital itself as well as the 
quadrupole moment, providing an evidence for the the- 
ory of field- induced multipoles in CeBg. We expect that 
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